DESY 08-206 
SFB-CPP/09-002 
January 2009 



Structural Relations of Harmonic Sums and 
Mellin Transforms at Weight w = 6 ^ 



Johannes Bliimlein 



Deutsches Elektronen-Synchrotron, DESY, 
Platanenallee 6, D- 15735 Zeuthen, Germany 



Abstract 

We derive the structural relations between nested harmonic sums and the corresponding 
Mellin transforms of Nielsen integrals and harmonic polylogarithms at weight w = 6. They 
emerge in the calculations of massless single-scale quantities in QED and QCD, such as 
anomalous dimensions and Wilson coefficients, to 3- and 4-loop order. We consider the 
set of the multiple harmonic sums at weight six without index {—1}. This restriction is 
sufficient for all known physical cases. The structural relations supplement the algebraic 
relations, due to the shuffle product between harmonic sums, studied earlier. The original 
amount of 486 possible harmonic sums contributing at weight w = 6 reduces to 99 sums 
with no index {— 1}. Algebraic and structural relations lead to a further reduction to 20 
basic functions. These functions supplement the set of 15 basic functions up to weight 
w = 5 derived formerly. We line out an algorithm to obtain the analytic representation of 
the basic sums in the complex plane. 



■•■Proceedings of the "Motives, Quantum Field Theory, and PseudodifFerential Operators", held at the Clay 
Mathematics Institute, Boston University, June 2-14, 2008 



1 Introduction 



Inclusive and semi-inclusive scattering cross sections in Quantum Field Theories as Quantum 
Electrodynamics (QED) and Quantum Chromodynamics (QCD) at higher loop order can be 
expressed in terms special classes of fundamental numbers and functions. Zero scale quantities, 
like the loop- expansion coefficients for renormalized couplings and masses in massless filed theo- 
ries, are given by special numbers, which are the multiple C-values [1,2] in the known orders. At 
higher orders and in the massive case other quantities more will contribute [3]. The next class of 
interest arc the single scale quantities to which the anomalous dimensions and Wilson coefficients 
do belong [4-6], likewise other hard scattering cross sections being differential in one variable 
z — L/ L given by the ratio of two Lorentz invariants with support 2; G [0, 1]. A natural way 
to study these quantities consists in representing them in Mellin-space performing the integral 
transform 

m{mN)= tdzz"" f{z) . (1.1) 

JO 

In the light-cone expansion [7] these quantities naturally emerge as moments for physical reasons 
with iV G N. Their mathematical representation is obtained in terms of nested harmonic 
sums [8-10] 

Si>,m = E ^-^Ss{k), So{k) = 1 , (1.2) 

k=l 

which form a unified language. This is the main reason to adopt this prescription also for other 
quantities of this kind. The harmonic sums lead to the multiple ^-values in the limit N ^ 00 
ioT b 1. In the latter case the harmonic sums diverge. ^ To obtain a representation which is 
as compact as possible we seek to find all relations between the harmonic sums. There are two 
classes of relations : 

i) the algebraic relations, cf. [11]. They are due to the index set of the harmonic sums only and 
result from their quasi-shuffie algebra [12]. 

ii) the structural relations. These relations depend on the other properties of the harmonic 
sums. One sub-class refers to relations being obtained considering harmonic sums at and 
integer multiples or fractions of A^, which leads to a continuation of A^ G Q. Harmonic sums 
can be represented in terms of Mellin-intcgrals of harmonic polylogarithms Hs{z) weighted by 
1/(1 ± ^) [13], which belong to the Poincare-iterated integrals [14]. The Mellin integrals are 
vahd ior N G R, A^ > A^o- Ftom these representations integration-by-parts relations can be 
derived. Furthermore, there is a large number of differentiation relations 

M[fiz)]iN)^M\\n^iz)fiz)\iN) . (1.3) 



We analyzed a wide class of physical single scale massless processes and those containing a 
single mass scale at two- and three loops [4-6] in the past, which led to the same set of basic 
harmonic sums and, related to it, basic Mellin transforms. Like in the case of zero scale quantities, 
this points to a unique representation, which is widely process independent and rather related 



^Due to the algebraic relations [11] of the harmonic sums one may show that this divergence is at most of 

0(ln'"(A^)), where m is the number of indices equal to one at the beginning of the index set. 
^Generalized polylogarithms and Z-sums were considered in [15]. 
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to the contributing Feynman integrals only. The representation in terms of harmonic sums is 
usually more compact than a corresponding representation by harmonic polylogarithms, since 
i) Mellin convolutions emerge as simple products; ii) harmonic polylogarithms are multiple 
integrals, which are usually not reduced to more compact analytic representations. The latter 
requires to solve (part of) these integrals analytically. In the case of harmonic sums the analytic 
continuation of their argument N to complex values has to be performed to apply them in physics 
problems. As lined out in Ref. [16-18] this is possible since harmonic sums can be represented 
in terms of factorial series [19] up to known algebraic terms. Harmonic sums turn out to be 
meromorphic functions with single poles at the non-positive integers. One may derive their 
asymptotic representation analytically and they obey recursion relations for complex arguments 
N. Due to this their unique representation is given in the complex plane. 

In the present paper we derive the structural relations of the weight w = 6 harmonic sums 
extending earlier work on the structural relations of harmonic sums up to weight w = 5 [18]. The 
paper is organized as follows. In Sections 2-6 we derive the structural relations of the harmonic 
sums of weight w = 6 of depth 2 to 6 for the harmonic sums not containing the index { — 1}. The 
restriction to this class of functions is valid in the massless case at least to three-loop order and 
the massive case to two-loop order. In Section 7 we summarize the set of basic functions chosen. 
The principal method to derive the analytic continuation of the harmonic sums to complex values 
of is outlined in Section 8 in an example. Section 9 contains the conclusions. Some useful 
integrals are summarized in the appendix. 



2 Twofold Sums 



The following w = 6 two-fold sums occur : S±5,i{N), ^±4,±2(iV), 5_3,3(^) and S3,3{N), 5_3 _3(iV). 
The latter sums are related to single harmonic sums through Euler's relation. 



(2.1) 



with aAb — sign(a) •sign(6)(|a| -|- |6|). For the former six sums we only consider the algebraically 
irreducible cases. In Ref. [18] the basic functions, which determine the harmonic sums without 
index { — 1} through their Mellin transform, up to w = 5 were found : 



w = 1 : l/(x-l) 
w = 2 : \n{l + x)/{x + 1) 
w = 3: Li2(a;)/(x± 1) 
w = 4: Li3(a;)/(a;+ 1), 
w = 5 : Li4{x)/{x ± 1), 
Ul{x)/{x±l), 



S,,2{X)/{X±1) 

Si,3{x)/{X±1), S2,2{X)/{X±1), 

[ln(a;)^i,2(-a;) - Ul{-x)/2]/{x ± 1) 



In the following we determine the corresponding basic functions for w = 6. 
In case of the double sums we show that they all can be related to 



M 



1 + x 



(iV) 



(2.2) 
(2.3) 
(2.4) 
(2.5) 

(2.6) 



(2.7) 
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up to derivatives of basic functions of lower degree and polynomials of known harmonic sums. 
The representation of S±5^i{N) read : 

'Li5(x)^ 



M 



x-1 



(N) - S,{N)C, + S2{N)C, - Ss{N)C3 + 3,(2 



ifM 



Lin fx) 



with 



and 



1 + x 

I dxg{x)[f{x)]+ = 
Jo 



(TV) + ^C5ln(2) - S6 - S.^{N)C, + S.2{N)C, 



dx[g{x) -g{l)]f{x) 



15, /•! , L 

se = — ln 2 C5+ / dz- 
Ib Jo 1 



Li5(z) 



+ z 



(2.8) 

(2.9) 
(2.10) 

(2.11) 



being one of the basic constants at weight w = 6. For the determination of the constants in 
the alternating case we use the tables associated to Ref. [10]. To express one of the sums given 
below we also give a second representation of S-5^i{N), 

Xi5(-x) - \n(x)U4(-x) +ln^(x)Li3(-x)/2" 



+ (_1)(a^+i)M 



x + 1 

- ln^(a;)Li2(-x)/6 - ln^(a;) ln(l + a;)/24 

x + 1 

-^C5 [S-i{N) - S,{N)] - |C| + la + f C5ln(2) - s. 



(AT) 



+(_l)(JV+i)M 



(TV) 



(2.12) 



The other two-fold sums are 

/ 4Li5(— x) — ln(x)Li4(— x) 



x-1 



(AT) 



(2.13) 



S-4,2(N) = (-l)^M 



4Li5(a;) — Li4(a;) ln(a;) 
1 + x 



(N) 



2SQ 3 1 S 

+2C3^-3(iV) - 3C4^-2(iV) + 4C5^-i(iV) + —CI - ^Ci - ^C5ln(2) + 4^6 



(2.14) 



\C2 [S-,{N) - S,{N)] - \CzS.3{N) + yC4^-2(A^) - ^C5^-i(iV) 



+ (_l)^+iM 



4Li5(— x) — ln(a;)Li4(— x) 



1 + x 



(TV) 



^4,2 (iV) 



.'^Cl + -a + K^H2)^As^ 
840^^ le^'^ -r S5 V ; b 

4Li5(a;) — ln(x)Li4(a;) \ 



-M 



x-1 



(2.15) 

(iV) + 2(3^3 (iV) - 3C4^2(iV) + 4C5^l(iV) 

(2.16) 
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M 



6Li5(-x) - 31n(a;)Li4(-x) +ln^(x)Li3(-a;)/2 
x-1 



(AT) 



Q 21 4^ 

--Cs [S.s{N) - Ss{N)] - -aS2{N) + -CsSiiN) 



1 



1 + x 



(2.17) 



+■ 



(TV) 



+ (-1) 



N+l 



M 



'3 ln(x) [^1,3(1 -x)- U] + In'(x) [^1,2(1 -x)- Cs] /2 



1 + x 



271 ,o 81,2 45, , ^ 
■^Cl + ^C| + ^C5ln(2)-6.a 



-^3,3^ = 3C4^2W-6C5^i(7V)-6M 



'S'i,4(l - a:) - Cs 
x-1 



(TV) 



M 



'ln2(a;) [5i,2(l - a;) - Ca] /2 + 31n(a;) [5i,3(l - a;) - C4] 



x-1 



5|(iV) + S,{N) 



In the above relations Nielsen integrals, [20], given by 
occur. The corresponding functions 5*1^^(1 — x) are given by 



(TV) 



(iV) 
(2.18) 



(2.19) 



(2.20) 



'^1,2(1 


-x) 


'^1,3(1 


— x) 


^1,4(1 


-x) 



-U^{x) + log(x)Li2(x) + - log(l - x) \og\x) + C(3) 

1 1 

-hu{x) + log(x)Li3(x) - - log^(a;)Li2(a;) - - log^(a;) log(l - x) + C(4) 

-Li5(a;) + ln(x)Li4(a:;) — ln^(a;)Li3(a;) + ^ ln^(a;)Li2(x) 



+-ln^(x)ln(l-x) + C5. 



(2.21) 



They are used to express the respective sums in terms of the Mellin transforms of basic functions 
and their derivatives w.r.t. A^. 

The algebraic relation for 5*3,3 (iV) can be used to express M[(Li5(x)/(a;— 1))+](A^). The Mellin 
transform in S'_3,_3(A^) allows to express >S'_4,_2(^) and S'_4,2(^) through (2.12). S4^2{N) and 
'S'-3,3(-^) do not contain new Mellin transforms. Therefore the only non-trivial Mellin transform 
needed to express the double sums at w = 6 is M[Li5(x)/(l + x)]{N). 

In some of the harmonic sums Mellin transforms of the type 



^hj-x) 
x±l 



(2.22) 



contribute. For odd values of A; = 2Z + 1 the harmonic sums Si^^(^k-i){N), S^(^k-i),i{^) and 
S-i^-i{N) allow to substitute the Mellin transforms of these functions in terms of Mellin trans- 
forms of basic functions and derivatives thereof. 
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For even values of k this argument applies to M[Life(— a;)/(l + x)]{N) but not to 
M[Life(— x)/(l + x)]{N). In the latter case one may use the relation 



1 Lifc(x^ 
2'=-2 1 - 



Uk{x) ^ Uk{x) ^ Likj-x) _^ Uk{x) 



1 — a; 



1 + x 



1 — X 



1 + x 



(2.23) 



Since in massless quantum field-theoretic calculations both denominators occur, one may apply 
this decomposition based on the first two cyclotomic polynomials, cf. [21], and the relation 
between Life(x^) and Lifc(±a;), [22]. The corresponding Mellin transforms also require half- 
integer arguments. In more general situations other cyclotomic polynomials might emerge. The 
relation 



2k-i 



M 



'Lifc(x^) 
x^-1 



+■ 



N 



M 



x-lj^_ 

/ Li.(-.r) 
x-1 

Jo 1 + x 



{N) + M 



'LiJx) 



+ M 



+■ 



(N) + M 



x+l J 



(TV) 



-X] 



x + l 



{N) 
(2.24) 



determines M[Life(-x)/(l + x)]{N). For k^2,A the last integral in (2.24) is given by 







1 U^jx^) 3 
dx- = C2ln(2)--C3 



l+x 

, Li4(a;2) 

dx 

1 + x 



f 

Jo 



2 2S 
-ln(2)C2' + 3(2(3-^(5. 



(2.25) 
(2.26) 



The corresponding relations for M[Li;i.(~^)/(l + 2^)](-^) ^'I's : 



M 



x+l 



(iV) = 



M 



Li4(— a;) 
x+l 



(TV) = 



-1m 

2 



+ M 



3^ 

+5& - 

-iM 



+ M 

1 



\x-l 

Li2(-a^) 
x-1 

^C2ln(2) 
|Li4(x)j 

Li4(-x) \ 
x-1 



+■ 



/N-1 
[ 2 



(iV) - M 



N 



(N) - M 



25 



+ M 

^h{x) 
x + l 

+ M 

114(0;) 
x + l 



Li2(a^) 
X — 1 

(N) 



Li4(a;) 



X 



(TV) 



(iV) 



(2.27) 



(TV) 



-^C2^1n(2)--C2C3 + ^Ca 



(2.28) 



In the case of w = 6 these relations do not lead to a further reduction of basic functions but are 
required at lower weights, cf. [18]. 
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3 Threefold Sums 



The triple sums are : 



'5-2,-3,1 (^) 
^l,-2,-3(iV) 

^l,-3,-2(iV) 
^-2,1,-3(A^) 



'5-3,2,1 (iV) 



-M 



x-1 J 



1 + X 



{N) + 5i(iV)(2C5 - C2C3) - jS2{N) + CsS,{N) 



(N) + (2C5 - C2C3)'5_i(iV) - ^S.2{N) 



(3.1) 



71 



29 



+C3'5-3(A^) + ^Ci + ^Ci - :^C,H^) - C2C3ln(2) + -s^ 



,5_3,-2,l(iV) = M 



840 
-^0,0,-1,0,1(3^ 
x-1 



32 



(iV) + C2-5-3,-i(iV) + [S.s{N) - Ss{N)] 



(3.2) 



C2ln(2)-|c3 



3 /Q fi7 

+^Cl^2(iV)-(-C2C3-Y^C5)5i(iV) 



(3.3) 
(3.4) 



S.2iN)S.,4N) + S,,,iN) + ^_3,-3(iV) - ^-3,1,-2 - ^-3,-2,1 
5„3(iV)^i,-2(iV) + ^i,5(iV) - S,{N)S_s,_2{N) - S_s,_s{N) 

+5_3,-2,i(iV) (3.5) 

S^{N)S_3,-2{N) + ^-4,-2(iV) + ^-3,-3(iV) - .5_3,l,-2(iV) - 5-3,-2,l(A^) (3.6) 
5_3,-3(iV) - 5_3(iV)5i,_2(iV) - ^i,5(iV) + 5i(iV)^_2,_3(iV) 

-S.2iN)S.s,,iN) - S,,i{N) + + Si{N)S.s,.2{N) 

+5_3,i,.2(iV) (3.7) 
" Ai{-x)/2 + Ss,2{-x)-S2,2i-x)ln{x)) 



M 



X — 1 



^C3-^C2ln(2) 



(AT) 

[S^s{N) - Ss{N)] 



-lciS2{N) 



23 7 ■ 

T7C5 — o^2C3 

.16 8 



(3.8) 



(-l)^M 



2Ss,2{x)-A^{x)/2 
x + l 



iN) + C2S-3,iiN)--aSMN) 



- (^C. - 3C2C3) S_.iN) + ^Cl + f,Ci + gCaln(2) + lc|ln^(2) 
+^C2C3ln(2) - ^C2ln^2) - 2C2Li4 (^) - 



(3.9) 



'52,-3,i(A^) = (-l)^M 



-f^o, -1,0,0, 1(2^) 



1 + x 

+C2S2MN) - C3S2,-l{N) 



(N) - ( ^C5 - f C2C3 ) S^,{N) 



+ 



-^Cl + 2Li4 Q + ^C3ln(2) - ^C2ln^(2) + ^ ln^(2) 



[S2{N) - S^2{N)] 



Si,2MN) = 

^l,-3,2(A^) = 



-^C2C3ln(2) - lc2ln^(2) - 4C2Li4 (i) + (2^111^(2) + 11(1 + 

-^C5ln(2)-^S6 (3.10) 

'5-3(A^)5i,2(iV) + Si,.5{N) - S,{N)S.s,2{N) - S.s,s(N) + S.s,2,i{N) (3.11) 

-S2{N)S_s,,{N) - -S_5,i(7V) + ,S2,_3,i(iV) + S,{N)S_s,2{N) 

+S-4,2{N) (3.12) 
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Si,2AN) 

S2,1,3{N) 
-51,3,2 (iV) 



S'3,-3(iV) - S.3{N)Si,2{N) - 5l,-5(iV) + Si{N)S2,-3{N) - S2,-3,l{N) 

+S2MN) + S,{N)S.sAN) + S.sAN) - ^-3,2,i(iV) (3.13) 
S2{N)S.3,,{N) + S_,,,{N) + S_s4N) - S2,-3,i{N) - 5-3,2,i(A^) (3.14) 



{3/2)A,{x)-U2{x)Us{x) 
x + 1 



{N) + C2S-2,2{N) - C3S-2,l{N) 



+ (^C5 - 2C2C3) S.,{N) - i^Cl - + |c5ln(2) - C,^ln^(2) 
+ ^C2C3ln(2) + ic2ln'(2) +4C2Li4 



-53,-2,1 (iV) 
^1,-2,3(A^) 



Sl,3,-2{N) = 

^-2,1,3(A^) = 
'53,l,-2(iV) = 



(3.15) 

,S3,_3(7V) - 53,i,-2(A^) - .5_2,3,i(A^) + S_2{N)Ss,,{N) + 5_5,i(A^) (3.16) 
>S_3,3(A^) - S3{N)S-2,i{N) - S-2,4{N) + S-2{N)Si,s{N) + S,,-,{N) 
+2S^2(N)Ss4N) + Ss,^s(N) - Si(N)Ss,-2(N) - S^M^) + S-5,i(N) 
-S3,i,-2{N) - S-2,3,l{N) (3.17) 
-S^2{N)S3,i{N) - S^5,i{N) + S-2,3,i{N) + Si{N)Ss,^2{N) 

+S4,-2{N) (3-18) 
S3{N)S^2,i{N) + S-2a{N) - S^2{N)Ss,,{N) - S3,-3{N) + -S3,i,-2(iV) (3.19) 



1)^M 



A,(-x)/2 + Ss,2{-x) - \n(x)S2,2(-x) 



1 + x 



(N) 



-^C2[S3,m-S3,.,{N)] 



^C3-^C2ln(2) 



[Ss{N) - S_s{N)] 
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lis 17 1 



-53,2,i(A^) = M 



'253,2(.x) - A^{x)/2' 



x — 1 
3 



(TV) 



H2S3,i{N) - \US2{N) - (H^5 - 3C2C3) S,{N) 
/ (3/2)^(x)-Li2(x)Li3(x) ^ 

V ^-1 J 



(3.20) 



(3.21) 



M 



(TV) 



-52,2,2(A^) = -M 



+C252,2(iV) - C3.52,l(iV) + (^^CS - 2C2C3) S^{N) 

53(iV)5i,2(iV) + 5i,5(iV) - 5i(iV)53,2(iV) - 53,3(iV) + ^3,2,l(iV) 
253,3(iV) - 53(iV)5i,2(iV) - 5i,5(iV) + Si{N)S2,3{N) - S2,3,l{N) 
+S2,a{N) + ^i(iV)^3,2(iV) - .53,2,l(iV) 

-^2(iV)^3,l(iV) - .55,l(iV) + 52,3,l(iV) + Si{N)S^,2{N) + ^4,2(iV) 
52(iV)53,i(iV) + ^5,l(iV) + .53,3(iV) - .52,3,l(iV) - ^3,2,l(iV) 

2Ai{x) + Ul{x) \n{x)/2 - 2^2,2(2;) ln(a;)^ 



(3.22) 
(3.23) 

(3.24) 
(3.25) 
(3.26) 



-M 



a; — 1 

'453,2(0;) -2Li2(a;)Li3(a;) 



(TV) 



x-l ) _ 

+2C352,l(iV)+2(C5-C2C3).Si(iV) 



(TV) 
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^2,-2,2(iV) 
^2,2,-2(A^) 



S. 



2,2,-2 



(AT) 



-52,-2,-2(A^) 
5-2,-2,2(iV) 



lsi{N) + ^S2iN)S,iN) + ^SeiN) 



(3.27) 



^-2,2,2(A^) = (-1)(^+^)M 
+(_l)(iv+i)M 



2Ai{x) + Ul{x) ln{x)/2 - 2^2,2(2;) ln(x) 
x+1 

^ASs,2{x)-2U2ix)Us{xy 



(TV) 



(iV) 



x + 1 

+2C3^-2,i(iV) + 2 (C5 - C2C3) ^-i(iV) 
-^C2C3ln(2) + ^C2ln^(2) + 2C2Li Q - ^CllA'^) 
23 

-— C5ln(2)+496 

-25_2,2,2(A^) + S2{N)S.2,2{N) + ^-4,2(A^) + S.2a{N) 
S-2,2,2m - \ [S2{N)S^2,2{N) + 5_4,2(A^) + ^-2,4(A^) 

-S2{N)S2,-2{N) - S,MN) - S2MN)] 

' -4.S3,2(-x) - \n(x)Ul(-x)/2 + 2 \n(x)S2,2(-x) ^ 

x-1 J 

'f 2Lh{-x)U2{-x) - 2Ai{-xy 

A ^-1 



37 
80' 



(3.28) 
(3.29) 

(3.30) 



M 



+ M 



(AT) 



(N) 



|C3^_2,-l(iV) + ic2 [^-2,-2(iV) - S_2,2m 



+ 



+ 2C3ln(2)-C2ln'(2) + -ln^(2) 



/II 23 

+ (xC2C3-YC5)^l(iV) 



[^-2(A^) - S2{N)] 

(3.31) 



- [->S_2,2,-2(A^) + S^2{N)S2,-2iN) + -5_4,-2(A^) + S24N)] 
^ [-5_2,2,-2(iV) + S.2{N)S.2,2{N) + S4,2{N) + ^-2,-4(A^)] 



(3.32) 
(3.33) 
(3.34) 



There emerge numerator functions, which do not belong to the class of Nielsen integrals ^, 

Aiix) = 
A2{x) = 



y 

^ln(l-y)5i,2(?/) 
y 

dy 
/o y 



A,{x) = r^[Li4(l-y)-C4] 

^0 V 



(3.35) 



As seen in Eqs. (3.27), {Ai{x)/{x — 1))+ is not a basic function since its Melhn transform 
reduces to single harmonic sums and known MelUn transforms algebraically. Furthermore, some 
numerator functions are given by harmonic polylogarithms i7ai,...,afe(a^), £ {~1)0, +1}, 
which cannot be reduced significantly further. Harmonic polylogarithms are Poincare-iterated 



*Note a misprint in Eq. (14), [17]. Li4(y) should read Li4(l — y). 
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integrals [14] over the alphabet [/o, /i, = [l/x, l/{x — 1), l/{x + 1)], [13], with 



Hi{x) 



ln(a;) 
— ln(l — x) 
ln(l + x) 



(3.36) 
(3.37) 
(3.38) 



and 



dy fa{y)H^{y) 



(3.39) 



4 Fourfold Sums 



The quadruple-index sums are : 

X + 1 



^-3,i,i,i(A^) = (-l)^M 



-52,-2,1,1 (iV) 



^-2,1,1,2 (A^) 



(TV) + CaS-2{N) - (2C5 - C2C^)S-r{N) 



1 1 / 1 \ 1 9^7 7 

+ gC2C3ln(2) - -C2ln^(2) - (2^4 (-) + -C|ln^(2) - ^Cl + ^C3ln^(2) 

" + 21n(2)Li5 (i) + ln^(2)Li4 

+lln^(2) + 2L.(i)-f 



-M 



<S'2,3(a;) 



x-l 



{N) + aS2{N)-{2C,-C2C3)Si{N) 



(4.1) 
(4.2) 



-1) 



Af+l 



M 



352,3(3;) + ^2(0;) 
x+1 



/II 



(iV) + C35-2,l(iV) + (^yCs - 3C2C3j ^-l(A^) 

41 1 7 
ggQS2 ^2^^^ ^ ' 



■^C2C3ln(2) + ic2ln^(2) + 2C2Li4 Q - ic|ln^(2) + 
|c| + ^C5ln(2)-41n(2)Li5 Q -21n^(2)Li4 Q - ^^A'^) 



'1\ 9 

-4Li6( - j + -S6 



(4.3) 



(-1) 
+ 



Af+l 



M 



^^0,-1,0,1,1 (a;) 



(iV) + (3^2,-1 (iV) + 



11 _ 41 ■ 
16 ^ - ^Cs^ 



5-1 (iV) 



-Li4 ( I) + la + ^C3ln(2) + lc2ln^(2) - ^ln^(2) [S,{N) - S.,{N)] 



-^C2C3ln(2) - ^C2ln^2) - 2(2^4 Q + ^(2' W(2) - ^Cl + ^C3ln^(2) 



' ^C5ln(2) +41n(2)Li5 (^) + 2 ln^(2)Li4 Q + ^^^\^) 



+4Li6 ( 2 j + «6 



(4.4) 



;-i)^M 



l + x 



{N) - (C2C3 - SCs) 5_i(7V) - (3^-2,1 (A^) + (2^-2,1,1 (A^) 



+^C2C3ln(2) + ^C2ln^2) + ^C2Li4 Q - IcIlA^) 



il 3 _ 27 , 
120^^ 16^^ 
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-2,-2,1,1 



(N) 



-3^C5ln(2) + -.e 

'-f^o,-i,o,i,i(.«) 



(4.5) 



-M 



X — 1 



(N) + (3^-2,-1 (A^) + (^C2C3 - ) S,{N) 



(4.6) 



^2,2,l,l(A^) 



M 



'3S2,3{x)+A2{xy 



X 



(N) + C3S2AN) + (yC5 - 3C2C3) S,{N) . 



Here, the harmonic polylogarithm i?o,-i,o,i,i(2^) is given by 

r-dy fy ,Si,2(z) 



rr t , rdy fy bi 

//o,-i,o,i,i(x) = J -I dz- 



+ z 



(4.7) 



(4.8) 



We tested the above sum-relations containing harmonic polylogarithms in the MeUin transforms 
numerically using the code of Ref . [23] . 

5 Fivefold Sums 



Two 5-fold sums contribute : 

'Si,,{xy 



'5-2,1,1,1,1 (iV) 



-M 



x-1 



1)^+1 M 



'5'i,4(a^) 
1 + x 



{N) + C,S,{N) 

{N) + (,S_,{N) 

1 . . . /1\ 1 



(5.1) 



+^C2C3ln(2) + ^C2ln^(2) + ^C2U, (^-j - -ain^(2) - ^C3ln^(2) 



49 



128 



a - ln(2)Li5 Q - \ \n\2)U, (^) - ^ ln^(2) - Lie (^) + ln(2)C5 ■ 

(5.2) 



All other sums can be traced back to these sums using algebraic relations [11]. The other Mellin 
transforms emerging in their representation were all calculated in Refs. [9, 18] before. 

6 Sixfold Sums 

Only one sixfold sum contributes at w = 6, -S'i^i^i^i^i^i(A^). This sum is completely reducible into 
a polynomial of single harmonic sums, cf. [9], 



^1 1 

6 



11 1111 

— -S,^ + - S2 St + Y^SsS! + - S,Sl + -S,S^ + - S^Sl 
+lSiS2S3 + j^Sl + ^S2S^ + ^ si + ^Se 



(6.1) 
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7 The Basic Functions 



In the following we summarize the basic functions the Mellin transforms of which represents the 
harmonic sums up to weight w = 6 without those carrying an index { — 1}. The corresponding 
sums of lower weight were determined in Refs. [8, 18,24]. The new 20 functions are given by 

w = 6: Li5(x)/(x + l), Si,4x)/{x±l), S2,3{x)/{x±l), 
Ss,2{x)/{x ± 1), U2{x)Us{x)/{x ± 1), 

Ai{x)/{x + l), A2{x)/{x±l), As{x)/{x+l) (7.1) 

^0 -i,o,i,i(^)/(^ ± 1)' ^0,0 -i,o,i(^)/(^ ± 1) 
[A,{-x) + 2S3,2{-x) - 2S2,2{-x) ln(x)]/(x ± 1) 

[Ai{-x) + 2S3,2{-x) - S2,2{-x) ln{x) + Ul{-x) ln(x)/4 - U3{-x)U2{-x)]/{x - 1) 

and extend the set Eqs. (2.2-2.6). The algebraic relations allow to express the initial set of 99 
functions by 30 functions and the structural relations reduce the basis further to 20 functions. 

8 Complex Analysis of Harmonic Sums 

The anomalous dimensions and Wilson coefficients expressed in Mellin space allow simple rep- 
resentations of the scale evolution of single-scale observablcs, which are given by ordinary differ- 
ential equations. The experimental measurement of the obscrvables requires the representation 
in ^-space. Therefore, one has to perform the analytic continuation of harmonic sums to com- 
plex values of N. Precise numerical representations for the analytic continuation of the basic 
functions up to weight w = 5 were derived in [25] based on the MINIMAX- method [26]. One 
may even obtain corresponding representations for quite general functions ^{z), z G [0, 1], as 
worked out for the heavy flavor Wilson coefficients to 2-loop order in [27]. ^ For other effective 
parameterizations see [28]. 

Here we aim on an exact representations. The inverse Mellin transforms are obtained by a 
contour integral around the singularities of the respective functions in the complex plane. 

We traced back all the harmonic sums to MeUin transforms of basic functions fi{z), 

Fr{N) = / dzMz) -, Ft{N) = / dzMz) ^ ^ . (8.1) 

Jo z — 1 Jo z + 1 

Eqs. (8.1) imply the recursion relations 

F-{N + 1) = -Fr{N)+ f'dzz'^Mz) , (8.2) 

Jo 

F+(N + 1) = F+(N) + (-1)^+1 dzz'^Mz) . (8.3) 

Jo 

The remaining integrals are simpler Mellin transforms, which correspond to harmonic sums of 
lower weight. 

^For another proposal for the analytic continuation of harmonic sums to € R, for which some simple 
examples were presented, cf. [29]. 
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If the functions fi{z)/{z — 1), fi{z)/{z + l) are analytic at 2; = 1 the MeUin transforms (8.1) 
can be represented in terms of factorial series [19]. Not all basic functions chosen above have 
this property. A corresponding analytic relation replacing 



Mz) ^ Ml - z) 



(8.4) 



always exists. The additional terms are lower weight functions in N or are related to these by 
differentiation. We use this representation and consider the factorial series. Due to this both the 
pole-structure and the asymptotic relation for \N\ — > 00 are known. ^. The poles are located 
at the integers below a fixed value Nq. The recursion relations (8.1) are used to express the 
respective harmonic sums at any value N E C except the poles. 

Let us illustrate this representation in an example for the harmonic sum <S'2,i,i,i,i(A^). The 
corresponding basic function is 



'Sia(z) 



(8.5) 



The recursion relation is given by 



M 



z-1 



(N + l) 



M 



'Si^ijz) 
z-1 



{N)+M[S,4z)] (N) 



with 



M [S,4z)] (TV) 



1 



(8.6) 



(8.7) 



cf. [33]. 

The numerator function possesses a branch-point at 2; = 1. The contributions related to 
terms ln'^(l — z)/{z ± 1) contained have to be subtracted exphcitely due to there logarithmic 
growth (to a power) for |A^| 00. This is either possible using the relation Si^4{z) to Li5(l — z) 

S,4z) = -Li5(l -z)+ ln(l - ^)Li4(l -z)-^ ln2(l - ^)Li3(l - ^) + ^ ^^H^ - z)U2{l - z) 



1 

24 



ln^(l-z)ln(z) + C5 



(8.8) 



or considering harmonic sums, which are algebraic equivalent to the above and are related to a 
basic function which is regular at 2; — > 1. We will follow the latter way and use the algebraic 
relations [11] to express 5'2,i,i,i,i(A'") afterwards, 



1,1 



'S'1,1,1,1,2 + ^ 



1 

12 

1 

'4 



'S'i5'2,i,i,i + jSa^i,!,!, + 'S'2,2,1,1 + 'S'2,1,2,1 + 'S'2,1,1,2 

>S'l>S'i^2,l,l + >S'2,2,l,l + 'S'1,3,1,1 + 'S'1,2,1,2 ~ 'S'i5'i^i^2,l ~ *S'2,l,2,l — >S'i^i^3_i — >S'i^i^2,2 



'S'l^S'i, 1,1,2 + 'S'2,1,1,2 + 'S'1,2,1,2 + 'S'l, 1,2,2 + >S'l, 1,1,3 



(8.9) 



''In [30] asymptotic relations for non-alternating harmonic sums to low orders in l/N^ were derived. Our 
algorithm given below is free of these restrictions. The main ideas were presented in January 2004 [31], see 
also [32]. 
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through known harmonic sums of lower weight. The latter sum obeys the representation 



Si,i,i,i,2{N) = 



M 



Lisfl - x) 



1 — X 

+C4-Si,i(iV)-C5^l(iV)+C6. 



{N)+C2Sl,^,,,^{N)-C^S,,^,,{N) 



(8.10) 



The function in the remaining Mellin transform is regular at ^ = 1 and can be represented in 
terms of a factorial series. The remainder terms in (8.10) are polynomials of single harmonic 
sums. Therefore the poles of 5'i,i,i,i,2(^), resp. 5'2,i,i,i,i(A^), are located at the non-positive 
integers. Finally we need the asymptotic representations of M [Li5(l — x)/{l — x)] (N), 



M 



Li5(l - x) 



x 



1 1 



179 1 



515 1 



216383 1 



183781 1 



7776 
4644828197 1 



+ 



41472 z* 
153375307 



194400000 z^ 25920000 z^ 
1 371224706507 1 



653456160000 z'^ ' 49787136000 z^ 25204737600000 z^ 
959290541 1 575134377343021 1 14855426650259 1 

^ + 



160030080000 z^'^ 16913534146740000 z^^ 312400053504000 
29106619674489691525729 1 225456132288901603 1 

319702820637227227200000 ^ ^ 



+ 



263567702701300558681 1 



788601079506240000 z^^ 
355061945309358701 1 



1053965342760089760000 z^^ 187184432058624000 z^^ 

1432477558547377054456843733 1 



4988266898917709221214400000 z^'^ 
192140702840923335916939 1 

Is 



13028192458306945920000 
2027981189268747465011536794768001 



1 



254294408120596135866406712880000 z^^ 



O 



y20 



(8.11) 



The corresponding representations for all other harmonic sums of weight w 
a forthcoming paper. 



6 will be given in 



9 Conclusions 

We derived the basic functions spanning the nested harmonic (alternating) sums up to weight 
w = 6 with no index { — 1}. This sub-class governs the functions contributing to the massless 
single-scale quantities, like the anomalous dimensions and Wilson coefficients to 3-loop order in 
QED and QCD. There arc first indications, that in the massive case, even in the limit ^ 
this class needs to be extended at 3-loop order, cf. [34]. Up to weight w = 5 all basic functions 
were given by polynomials of Nielsen integrals, Eq. (2.20), of argument x or —x weighted by 
l/(x ± 1). Although most of the basic functions at w = 6 share this property, some contain 
1-dimensional integrals over polynomials of Nielsen integrals 74j(±x)|^ 3 and more dimensional 
integrals, which are not reducible. This is generally expected and the cases up to w = 5 form an 
exception. 

We lined out how the exact the representation of the Mellin transforms of the basic functions 
can be obtained, generalizing effective numerical high-precision representations [25,27]. Up 
to terms which can be determined algebraically the Mellin transforms of the basic functions 
are factorial series. The singularities of the Mellin transforms are located at the non-positive 
integers. They obey recursion relations for N + 1. The asymptotic representation of 
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the Mellin transforms can be determined analytically. The basic Mellin transforms are thus 
generalizations of Euler's ■^-function and their derivatives, which describe the single harmonic 
sums. 
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10 Appendix A: Useful Integrals 

In this appendix we list useful constants and integrals. 

Lifc(l) = Ck (10.1) 

s,4i) = Cfc+i (10.2) 

-^2,2(1) = ^Cl (10.3) 

5-3,2(1) = 2C5-C2C3 (10.4) 
29 1 

Ss,2(-1) = -^C5+2C2C3 (10.5) 

52,3(1) = 2C5-C2C3 (10.6) 

A,il) ^ -3C5 + 2C2C3 (10.7) 

^l(-l) = -^C5 + ^C2C3 (10.8) 

^(1) = —Cs (10.9) 

^(1) = -3C5 + C2C3 (10.10) 

rdy^^^^tyl ^ ln{l + x)Us{-x) + lul{-x) (10.11) 
Jo 1 + y 2 

r^y Hy)U,(-y) ^ i^i + ^)i^^)Li2(-x) + ^Li^(-x) 
Jo i. + y 2 

-2S2,2{-x) +2\n{x)Si,2{-x) (10.12) 

r^y^lM. ^ ln{l-x)Si,2{x)+3Si,s{x) (10.13) 
Jo y-l 

r^[Li2(l-|/)-C2] = -2Li3(a;)+ln(a;)Li2(x) (10.14) 
Jo y 

r^y^Mvl ^ lLi2(a;)+ln(l-a;)Li3(a;) (10.15) 

Jo 1/ — 1 2 

rdy^-^U2{y) = ^Li^(x)+ln(x)ln(l-x)Li2(x) 

Jo y — I 2 

-2S2,2{x) + 2 \n{x)Si,2{x) (10.16) 

r^\n{l-y)Us{y) = -U2{x)Us{x) + Ai{x) (10.17) 
Jo y 

r^uMHy)H^-y) = -\h^)^^1{^) + \m^) (lo.is) 

Jo y 2 2 

r^Li2(-y)ln(y)ln(l + y) = -Un{x)Ul{-x) + ]-A,{-x) (10.19) 
Jo y 2 2 

r^Us{-y)\n{l + y) = -U2{-x)Us{-x) + A,{-x) (10.20) 
Jo y 

r^[Lk{l-x)-C3] = -S,,2{x)ln{x) + lul{x)-C2U2{x) (10.21) 

JO iJ ^ 
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r^S.MHy) = S2,2{x)]n{x)-Ss,2{x) (10.22) 
^0 y 

r^S,,2{-y)Hy) = S2,2{-x)ln{x)-Ss,2{-x) (10.23) 
JO y 
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